Spatial and polarization structure in 
g: micro-dome resonators: effects of a 

Bragg mirror 



(N 

3 



(N 

CZ3 



David H. Foster and Jens U. Nockel 
Oregon Center for Optics, 

1274 University of Oregon, 
gn; Eugene, OR 97403-1274 

g ■ http : //darkwing . uoregon . edu/~noeckel 

•i-H 

>v Published in Laser Resonators and Beam Contol VII, edited 

by Alexis V. Kudryashov and Alan H. Paxton, Proceedings of 

SPIE 5333, 195-203 (2004) 
^ ■ 

m 

\Q . Micro-domes based on a combination of metallic and dielectric multilayer 

mirrors are studied using a fully vectorial numerical basis-expansion method 
that accurately accounts for the effects of an arbitrary Bragg stack and can 
efficiently cover a large range of dome shapes and sizes. Results are examined 
from three different viewpoints: (i) the ray-optics limit, (ii) the (semi-) con- 
focal limit for which exact wave solutions are known, and (iii) the paraxial 

^ approximation using vectorial Gaussian beams. 
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1 Introduction 

Microresonators based on planar distributed Bragg reflectors (DBR) are ubiquitous be- 
cause they afford great design flexibility, e.g., tailored stopbands of ultrahigh reflectivity, 
and at the same their fabrication is well-developed in many material systems. In many 
applications, however, control over transverse mode profiles and sidemodes is improved 
when non-planar structures are used. Stable, dome-shaped cavities have been employed 
with InGaAs quantum- well lasers [1], and as passive filter cavites [2]. Such cavities are 
expected to be of great promise for semiconductor-based quantum optics [3], because 



strong focusing and large modulation of the local density of states can be achieved 
[4] . Because dome shaped DBR structures are technologically challenging, combinations 
of flat DBR stacks and curved metal mirrors have been considered as a compromise 
solutionfl]. Planar metal layers have recently proven useful as top mirrors in thin-film 
organic lasers [5] , and novel fabrication schemes [6] are likely to make curved metal mir- 
rors a building block in microcavity design. In this paper, we address dome cavities in 
which metal mirrors and Bragg stacks are combined. 

Computationally, such cavities are nontrivial even when the medium is linear, because 
they combine an unconventional three-dimensional shape with boundary conditions that 
differ at the top and bottom of the structure. Realistic Bragg mirrors of finite depth 
may also permit significant leakage, which has to be modeled accurately in order to 
determine the coupling properties of the resonator. We have performed fully vectorial 
electromagnetic simulations of axially symmetric cavities formed by a curved metallic 
dome on a stratified Bragg stack. The composition of the stack mirror and the shape 
of the metal dome have been varied in a wide range of parameters, and the capabilities 
of the numerical methods have been tested up to cavity lengths of 50A, where A is the 
wavelength in the dielectric comprising the bulk of the dome. 

2 Cavity geometry 

The cavity geometry is shown in Fig. 1. It is useful to define a cylindrical coordinate 
system (p, 0, z) with z being the symmetry axis. The height h of the dome above the 
Bragg mirror, and its (not necessecarily constant) radius of curvature, R, determine 
the type of modes the cavity can support. For definiteness, in this work the bulk of 
the dome and the space outside of the structure have refractive index no = 1. In the 
numerical results shown below, we will adjust the the Bragg stack design such that its 
stop band is centered at the modes of interest. The DBR layer structure begins with 
a spacer layer of optical thickness X s and refractive index 3.5, followed by N s = 20 
quarter-wave layer pairs AB of refractive index n^ = 3.0, n# = 3.5. Our main goal is to 
demonstrate the importance of the polarization-dependent phase shifts induced by this 
mirror in determining the resonator modes. Therefore, we also performed comparative 
calculations with the Bragg stack replaced by a conducting surface. In accordance with 
this main goal of the present study, absorption in the metal surfaces was neglected. 
Below, the curved mirror will be taken to be a spherical shell of constant R. However, 
we have also carried out calculations for domes in the shape of a rotational paraboloid, 
to compare with earlier calculations based on a separation of variables that is possible 
for that system [ i] when only scalar fields are considered. At size parameters between 
kR = 26 ... 31, and using only metal mirrors, all modes agreed in their spectral positions 
to better than five decimal places, despite the very different computational approach used 
in the present work. 
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Figure 1: Geometry of the dome resonator with a Bragg mirror at its base. The dome 
mirror is taken to be a spherical shell of radius R, centered at a depth R — h 
below the planar mirror surface. The bulk of the dome, and the substrate have 
refractive index n®. 

3 Vector fields in the dome resonator 



We are looking for the divergence- free, monochromatic (ex exp[— tut}) electric fields 
satisfying the vectorial wave equation. Using the axial symmetry, the field components 
E p , E^ and E z can be assumed to have a ^-dependence of the form exp(zm0). For fixed 
azimuthal mode number m, the electric field transverse to the z axis, E^ = E p p - 
can be decomposed into circularly polarized components using 



where 



With this, one has 
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The standard conducting-mirror boundary conditions on the curved dome (E\\ = 0, 
H± = 0), together with the planar mirror at the base of the dome, generally couple E p , 
E^ and E z . All modes can be labeled by m, and the substitution m — > —m (for m^O) 
leads to a degenerate mode. In presenting results later on, there is no loss of generality 
if m > is assumed. 

As was found in [1], focusing on the z axis is most pronounced for the states of the 
cavity in which either the magnetic or the electric field is polarized exlusively in the 
azimuthal direction 0. If B oc 0, then E has components only along z and p. 



4 All-metal cavity 

In a recent free-space beam experiment [7], extremely strong focusing was indeed ob- 
served with light having this type of polarization; the beam was called "radially polar- 
ized", but a significant contribution to the focused spot actually comes from E z . The 
analogue of this for our dome cavity is illustrated in Fig. 2, assuming the flat mirror to 
be a conductor 1 . In the hemispherical limit, this mode evolves smoothly into an even 
more strongly focused electric multipole with m = and 1 = 1. 
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Figure 2: Grayscale plots of the E-field components (a) E x and (b) ^ for a "radially 
polarized" mode at A = 867nm in a conducting dome with h = 9.95 /im, 
i? = 10 /im. Tight focusing is apparent, especially in E z . Here and in the 
following vertical cross-sections, the field is plotted in the xz plane. The 
projection of the electric polarization in a transverse (xy) plane at height 
z = 0.25 /im above the planar metal mirror are shown in (c). 



In the generic cavity with h < R, any given mode will have admixtures of many 
multipoles with different I, both of electric and magnetic type. The azimuthal mode 
number m is a fixed parameter in the calculation. To obtain Fig. 2, the coefficients y 
of this basis expansion are determined from a linear system of equations set up by the 
boundary conditions. The number of unknowns is determined by the number of different 
I — Imin ■ ■ • ^maxi n the expansion. Here, l m in — max(l, \m\), and an approximate cutoff 
Imax is given by the semiclassical limit knop max where p max is the maximum radial extent 
of the dome. Since we are interested in cavity modes, there is no incoming wave outside 



1 In all wave plots that are shown, we set u)t = it/2 or and take the real parts of the electric fields. 



the dome, and the resulting system of equations is initially homogeneous. Our method 
proceeds by adding a further "seed" equation that sets a certain linear combinations of 
unknowns equal to one [8], thus converting the problem into an inhomogeneous matrix 
equation Ay = b where the matrix A depends on k, and b ^ because of the seed 
condition. The boundary conditions are enforced by point matching (real-space points 
on the curved mirror, and k-space points for the planar mirror), and the wavenumbers 
of the cavity modes are found by minimizing the residual yply — b| 2 as a function of 
k. Details of the numerical method are presented in [9]. Before discussing the additonal 
complications posed by a dielectric stack mirror, we use the all-metal cavity to further 
investigate the effect of the dome-shaped mirror. 

5 Relation to ray dynamics and confocal limit 

A global picture of the possible cavity modes can be obtained from geometric optics. 
Mirror configurations with h < R are called stable because they support modes cen- 
tered on stable axial rays, while h = R creates a confocal geometry if one unfolds the 
cavity about the planar mirror [ 1] (one might alternatively refer to h = R as the "semi- 
confocal" condition, but ray-optically there is no difference to a symmetric, confocal 
mirror arrangement). This marginally stable limiting case features highly focused but 
non-paraxial modes. Non-paraxial ray orbits that can support cavity modes are also 
found for h < R, and an efficient tool to reveal all the coexisting types of stable mode 
patterns is the Poincare surface of section [4]. 

The axial symmetry of the cavity implies that any given ray trajectory r(s), where s 
is the path length, can be parametrized by an equation of the form r x dr/ds = L with 
a vector L that changes upon reflection at the boundaries but whose component L z is 
a constant of the motion (the skewness of the rays). Thus L z labels families of rays, and 
fundamental Gaussian-beam type modes are built on rays with L z m 0. However, even 
when L z is fixed, stable periodic ray orbits with a multitude of different topologies will 
generally be found. 

Regarding the internal ray dynamics, the dome cavity is thus closely related to axially 
symmetric spheroidal resonators, such as droplets and microspheres [10], for which the 
confocal limit corresponds to perfect sphericity. This analogy extends to the occurence 
of whispering-gallery type patterns and their non-perturbative breakup under shape 
distortion, which has recenty been observed in fused-silica microspheres [11]. As an 
illustration of this phenomenology in a dome cavity, Fig. 3 shows a cavity with h = 9 (im, 
R = 10 fim. The rectangular ray pattern in (a) for L z = can be thought of as a 
stabilized periodic whispering-gallery trajectory. 

A corresponding solution of the vectorial wave equation is shown in Fig. 3 (b), in- 
dicating that the potential for strongly focused modes in near-confocal cavities comes 
at the price of an increased variety of off-axis modes which need to be characterized in 
order to choose the optimal design. It is conceivable that such new mode patterns are 
in fact desirable for certain applications. In order to take advantage of this new variety 
of modes, it will be important to combine 3D cavity shape design with suitably tailored 




Figure 3: (a) ray pattern (b) y-component of the electric field (grayscale), plotted in a 
meridional cross section along the xz plane. Axis labels are in microns. 

Bragg mirrors that discriminate against the unwanted type of modes. For paraboloidal 
domes, some initial exploration of this aspect has been performed [4], in particular ap- 
plying ray-optic phase space analysis to the problem of angle- dependent Bragg mirror 
transmission. 



6 Relation to the paraxial approximation 

However in the remainder of this paper, we turn our attention to the more conventional 
modes centered on the z axis. These can be understood in terms of Gaussian beams 
when the paraxial parameter p = X/zr <C 1, where zr = ttwq is the the Rayleigh range 
and Wo is the beam waist radius [12]. An expression for the transverse field at a fixed 
m in terms of the vector Laguerre-Gaussian basis can be written in analogy to (2), 
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Here, A and B are amplitude coefficients, the scalar Laguerre- Gauss functions LG 7 ^ * 
have a dependence exp(?(m ± l)</>), and N is the order of the Laguerre- Gauss beam, 
N = \m\ + 1, 1 77i | +3, .... This definiton is related to a convenient normalized expression 
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fixed m). The special choice N = \m\ — 1 is also possible, provided m ^ 0; in this case, 
admissible solutions must have 74 = 0ifm>0, orl? = 0ifm<0, so that E^ has 
circular polarization. As a consequence, free-space LG solutions at given m with order 
N = 1 777 1 — 1 are nondegenerate. For all other orders N (and any 777), the two terms in 
(3) are degenerate free-space solutions in the paraxial approximation, and will remain 
nearly degenerate in a paraxial cavity. The polarization properties of the mode are fixed 
by the way in which the boundary conditions split this doublet (for 777 7^ 0). 

The fundamental Gaussian modes of the dome can be obtained by setting m — 1, N — 
and A = 0. The spatial dependence is then E? oc LGq, which is also called a TEMoo 
mode. Figures 4 (a) and (b) show how the resemblance to this Gaussian mode emerges 
starting from the hemispherical limit. 



Raising N in increments of 2 makes higher transverse excitations of the Gaussian mode 
accessible. This is again borne out by the numerical vectorial solutions, as shown for 
m — 1 in Fig. 4 (c) and (d). These modes are well described by (3) with N — \m\ + 1 — 2 
and either A = (c) or B = (d). Note that the all-metallic boundary conditions have 
uniquely fixed the polarization of the modes in Fig. 4 to be circular, i.e. only one of the 
terms in (3) is present in all cases shown up to this point, and the degeneracy in A is 
lifted by a small but finite amount. 




Figure 4: Some modes of the all-metal cavity with azimuthal order m = 1: fixing R = 
10 /iin, a mode of type TEMoo narrows as h decreases from (a) h = 9.9992 
(A = 889.02 nm) to (b) h = 9.9000 (A = 881.30 nm). Compare to the analogous 
results with a Bragg stack in Fig. 5. Shown in (c) is a paraxial mode of type 
LG\ at h = 10.5 /im, R = 70 /im and A = 796.060nm. (d) shows its near- 
degenerate partner LG\ at A = 796.055nm. See Fig. 6 for the transverse 
polarization patterns corresponding to (c) and (d). 



7 Combining metal dome and Bragg mirror 

The previous results raise the question if all cavity modes of the axially symmetric dome 
necessarily factor into a unique polarization vector multiplied by a spatial wave function, 
thus making the wave problem essentially scalar. We will demonstrate below that this 
is not the case when a Bragg stack replaces the planar metal mirror at the base of the 
dome. To this end we first discuss the necessary extension of the numerical method 
based on the basis expansion of Section 4, allowing us to treat the combined problem of 
metal dome and dielectric mirror without approximations. 

The complex reflection amplitude of an arbitrary planar Bragg grating is calculated 
with the transfer matrix method [14] and can be written in the plane-wave basis as a 
function of polarization, wavelength and angle of incidence, T sp (/c, 9). Here, subscripts s 



and p distinguish between polarization perpendicular and in the plane of incidence. Be- 
cause of this simplicity, one may decide to solve the whole dome problem in a discretized 
plane wave basis, which after symmetrization according to the axial rotation invariance 
leads to expansions of all vector field components in Bessel beams. Their (scalar) form 
in the dome region is 

ipBessei(p,4>i z i@) — 27rz m exp(zm0) exp (in k z cos 9) J m (pn k sin. 6), (4) 

where 9 is a parameter that specifies the cone of plane waves from which this beam is 
constructed. 

However, the approach we follow in this paper is to combine Bessel waves and mul- 
tipoles. We use the 4>Bessei to describe the DBR mirror, because reflection amplitudes 
for a Bessel beam of given 9 are identical to those of a plane wave at that same incident 
angle, r s>p (9). But we retain multipoles as the basis in which to specify the dome bound- 
ary conditions, because these basis functions evolve into the true cavity modes in the 
limiting case of the all-metal hemisphere. We discuss elsewhere [9] the relative merits of 
this approach compared to the procedurally more direct method of using a single basis 
for the entire domain. Our two-basis method makes it necessary to (i) discretize the 
cone angles 9 used in the expansions, and (ii) to implement a transformation between 
multipoles and plane waves so that we effectively obtain a formulation of the Bragg 
reflection in the multipole basis. 

With these modifications, the solution algorithm is the same as described for the all- 
metal cavity in Section 4. An important additional piece of information that enters for 
the Bragg-bounded dome is the Q- factor of the modes: if the Bragg reflectivities r s>p {9) 
are not of modulus one, there is leakage loss which makes the modes metastable; to find 
these modes, their wavenumbers k are now allowed to be complex [15]. In the following 
we consider only modes calculated to have Q > 10 4 , but omit a more detailed discussion 
of the Q-factors because that will require addressing the competing absorption losses in 
the curved metal mirror as well. 

Figure 5 lists a set of results analogous to the ones shown in Fig. 4, but now including 
the effect of the Bragg mirror. There is no significant distinction between the two systems 
from a purely ray-optic point of view, and the modes in parts (a,b) of Figures 4 and 
5 were obtained in the same way, by adiabatically following a fundamental Gaussian 
mode from small to large h. Nevertheless, the near-hemispherical mode of Fig. 5 (a) 
exhibits a V-shaped distribution in E y which is absent for all-metal mirrors. We have 
found the same behavior in cavities as large as R = 40 /mi, h = 39.7 /mi. The mode is 
predominantly s-polarized, and its formation can be interpreted as a consequence of the 
fact that the Bragg reflectivity r s (9) has a phase that depends on angle of incidence; in 
this sense, this phenomenon is induced by the Bragg stack, but arises essentially as a 
scalar wave effect [9]. 

The vector nature of the electromagnetic field enters in a more intricate way if we 
compare the modes shown in Fig. 5 (c), (d) to the corresponding ones in Fig. 4 (c), (d). 
For this cavity shape, all modes shown are well approximated by Eq. (3) with m = 1 
and N = 2. However, whereas A = or B = in Fig. 4, we now find A, B ^ in 




Figure 5: Modes with m = 1 in the cavity bounded by metal dome and dielectric Bragg 
stack. The cross sections in (a) and (b) were obtained with R = 10 /mi for 
modes at dome heights (a) h = 9.782 /mi and (b) h = 9.544 /mi. The calculated 
wavelengths and Q-factors are (a) A = 812.394nm, and (b) A = 796.328nm. 
A striking difference between (a) and the comparable metallic-mirror mode 
in Fig. 4 (a) is the V-shaped spatial pattern; this disappears at the slightly 
more paraxial h in (b). (c) and (d) show Gauss-Hermite like modes for the 
same cavity shape as in Fig. 4 (c), (d). The modes form a doublet at (c) 
A = 811.063nm and (d) A = 811.068nm. 



both members of the doublet. Since the spatial dependence of LG\ and LG\ is different, 
polarization and orbital part of the wave are entangled in the sense that it is impossible 
to factor out a uniform polarization vector. 

The difference between the side views in Fig. 5 (c), (d) and Fig. 4 (c), (d) is due 
to the fact that the spatial mode profiles in the dome with Bragg stack more closely 
resemble Hermite-Gauss beams[l3], u^ . This can be reconciled with Eq. (3) by noting 
that both LG\and LG\ have projections onto Uq§ an d i*fo ■ 

However, a transverse cross section of the mode field reveals the entanglement of 
polarization and spatial structure in these modes, which is not obvious from the side 
views alone. In Fig. 6, we compare the paraxial modes from Figs. 5 and 4 with respect 
to the instantaneous projections of the mode field into a horizontal plane. Because Eq. 
(3) is approximately valid, one can in fact construct [9] both (c) and (d) by appropriately 
superimposing the circularly polarized fields (a) and (b). Neglecting the z components 
of the electric field, the resulting modes in Fig. 6 (c) and (d) are close to circular 
polarization near the z axis, but generally exhibit varying polarization over the cross 
section. The vortex pair bracketing a central high-field region in Fig. 6 (c) corresponds 
to a ring of radial (linear) polarization when the time dependence of the wave is taken 
into account. 

It is remarkable in the Bragg-mirror dome that the particular transverse polarization 




Figure 6: Transverse polarization cross sections, (a) and (b) correspond to the modes 
in Fig. 4 (c), (d); the section is taken at height z = 0.25 /im and is cropped 
to contain only the regions of significant field strength. To see that this repre- 
sents circular polarization, imagine the pattern rotating rigidly, and track the 
resulting arrow direction at any fixed point. Similarly, (c) and (d) correspond 
to the modes in Fig. 5 (c), (d), but plotted at z = 0.05 /mi to intersect the 
domains of strongest field. 

patterns in Fig. 6 (c), (d) are in fact robust over a large range of cavity shapes h and 
also quite insensitive to the spectral location of the mode within the stop band of the 
Bragg mirror. 

The Laguerre-Gauss modes responsible for the persistent mixing with polarization 
structure given by Fig. 6 (c), (d) have order N = 2. For completeness, we verify that 
the other remaining mode of order 2, which has m — 3, does not undergo the same 
qualitative transition when changing from metal to Bragg mirror. The reason is that in 
this case N = \m\ — 1 and thus there is no doublet from which superpositions can be 
formed. The circularly polarized singlet modes of the paraxial dome remain circularly 
polarized independently of the nature of the planar mirror. This is shown in Fig. 7. 
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Figure 7: Transverse polarization cross sections at height 2 = 0.25 /im for a mode cor- 
responding to m = 3, iV = 2. This is a circularly polarized singlet state 
because N — \m\ — 1. Consequently, the polarization patterns with metal mir- 
ror (a) at A = 796.091nm and Bragg stack (b) at A = 811.086 are practically 
indistinuishable. Insets show the corresponding profile of E x in the yz plane. 

8 Conclusion 

There exists a great wealth of other stable but non-paraxial modes in the dome cavity 
with h < R, which can be classified with the help of ray-based methods as exemplified by 
Fig. 2. The exact numerical calculations demonstrated in the present work are relevant 
to such quasiclassical studies because we are able to treat realistic cavities whose size is 
large in relation to A. Future work in this direction will in particular have to address 
complex boundary conditions as they arise when different types of mirrors make up the 
3D cavity. In our case, a combination of conducting surface and dielectric multilayers is 
considered. 

The vectorial nature of the cavity fields is essential in this system. Fully vectorial 
mode calculations of the type performed here indicate that one of the effects induced by 
a Bragg-stack is the persistent mixing of doublets illustrated in Fig. 6. The nonuniform 
polarization patterns can be further analyzed, in particular regarding the locations of 
their singularities: in the near-paraxial situation, we pointed out the occurence of linear 
polarization on rings surrounding the circularly polarized beam axis. These considera- 
tions are of significance in particular when the coupling between cavity field and dipole 
emitters at the base of the dome is considered [4] . 

Although we presented only results for dome mirrors in the form of a conducting 
spherical shell, other shapes such as paraboloids can be treated. By making contact with 
various limiting cases in this paper, the general numerical techniques have been validated. 
A more detailed discussion of the numerical methods and their implementation can be 
found in a forthcoming publication [9]. 



11 



ACKNOWLEDGMENTS 

This work was supported by NSF Grant ECS-02-39332. 

References 

[1] F. M. Matinaga, A. Karlsson, S. Machida, Y. Yamamoto, T. Suzuki, Y. Kadota, and 
M. lkeda, "Low-threshold operation of emispherical microcavity single-quantum- 
well lasers at 4 k," Appl. Phys. Lett. 62, pp. 443-445, 1993. 

[2] M. Aziz, J. Pfeiffer, and P. Meissner, "Modal behaviour of passive, stable micro- 
cavities," Phys. Stat. Sol. (a) 188, pp. 979-982, 2001. 

[3] T. M.Stace, G. J. Milburn, and C. H. W. Barnes, "Entangled two-photon source 
using biexciton emission of an asymmetric quantum dot in a cavity," Phys. Rev. B 
67, p. 085317, 2003. 

[4] J. U. Nockel, G. Bourdon, E. L. Ru, R. Adams, J.-M. M. I. Robert, and I. Abram, 
"Mode structure and ray dynamics of a parabolic dome microcavity," Phys. Rev. E 
62, pp. 8677-8699, 2000. 

[5] V. Bulovic, V. G. Kozlov, V. B. Khalfin, and S. R. Forrest, "Transform-limited, 
narrow-linewidth lasing action in organic semiconductor microcavities," Science 
279, pp. 553-555, 1998. 

[6] S. Coyle, G. V. Prakash, J. J. Baumberg, M. Abdelsalem, and P. N. Bartlett, 
"Spherical micromirrors from templated self-assembly: Polarization rotation on the 
micron scale," Appl. Phys. Lett. 83, pp. 767-769, 2003. 

[7] R. Dorn, S. Quabis, and G. Leuchs, "Smaller, sharper focus for a radially polarized 
light beam," Phys. Rev. Lett. 91, p. 233901, 2003. 

[8] J. U. Nockel, PhD thesis, Yale University, 1997. 

[9] D. H. Foster and J. U. Nockel, "Methods for 3-d vector microcavity problems in- 
volving a planar dielectric mirror," Opt. Commun. 234, pp. 351-383, 2004. 

[10] A. Mekis, J. U. Nockel, G. Chen, A. D. Stone, and R. K. Chang, "Ray chaos and 
q-spoiling in lasing droplets," Phys. Rev. Lett. 75, pp. 2682-2685, 1995. 

[11] S. Lacey, H. Wang, D. Foster, and J. Nockel, "Directional tunnel escape from nearly 
spherical optical resonators," Phys. Rev. Lett. 91, p. 033902, 2003. 

[12] A. E. Siegman, Lasers, University Science Books, Sausalito, CA, 1986. 

[13] M. W. Beijersbergen, L. Allen, H. E. L. O. van der Veen, and J. P. Woerdman, 
"Astigmatic laser mode converters and transfer of orbital angular momentum," 
Opt. Commun. 96, pp. 123-132, 1993. 



12 



[14] P. Yeh, Optical waves in layered media, Wiley, New York, 1988. 

[15] J. U. Nockel and R. K. Chang, "2-d microcavities: Theory and experiments," in 
Cavity-Enhanced Spectroscopies, R. D. van Zee and J. P. Looney, eds., Experimental 
Methods in the Physical Sciences 40, pp. 185-226, Academic Press, San Diego, 2002. 



13 



